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1. Introduction 

Einstein General Relativity is accepted as a correct theory of gravity at huge range of 
scales, from cosmological to the subatomic ones, that successfully describes almost all 
classical gravitational phenomena (with few exceptions like singularities, dark energy, 
etc) . In spite of this, there is no consistent theory of quantum gravitational phenomena, 
that is of physics at very small length scales. It is expected that the general relativistic 
description of gravity, and, as the result, of the space-time, is inadequate at short 
distances (and maybe at very large distances due to the cosmological constant problem). 
There are many different proposals how to extend General Relativity. Aside from the 
well established approaches such as string theory and loop gravity some work has been 
done in noncommutative theory of gravity associated with noncommutative spaces 
H, I, I, i as well as "matrix gravity" g, |, |3, |, |, |10|, [U], |12|, |l| . 



One of the first attempts of constructing matrix gravity was proposed by Okubo 
and Maluf 0, whose primary goal was to unify Einstein's general relativity and Yang- 
Mills gauge theory. The main ingredient in their approach was the matrix valued affine 
connection and the main purpose of their work was the unification of Einstein general 
relativity and Yang-Mills theory. They introduced the matrix-valued metric (or the 
frame vectors) and the connection coefficients and developed the formalism of some 
kind of "matrix-valued" differential geometry, that is, covariant derivatives, curvatures 
etc. Because of the non-commutativity many choices have to be made in any "matrix" 
differential geometry, such as how to raise indices, how to define covariant derivatives, 
what constraints to use to fix the torsion and some others. Different choices lead 
to different physical theories and there is no physical principle which could help in 
distinguishing a unique possibility. The problems encountered by Okubo and Maluf 
were: the presence of the spin-3 components of the connection, non-uniqueness of con- 
nection, non-covariance of the torsion, non-covariance of higher covariant derivatives, 
consistency of covariant derivatives of higher rank tensors and some others. The further 
work in this direction [B| was focused on an particular ansatz for the affine connection 



which was introduced originally by Einstein and Kaufmann [|14[ in the Abelian theory. 
This ansatz separates the usual Christoffel coefficients and the Yang-Mills gauge fields 
and leads to a standard Einstein general relativity coupled to the Yang-Mills theory 
with a particular gauge group. There are no new physical degrees of freedom in this 
formulation. 

The formalism of matrix geometry is related to the algebra-valued tangent space 
formulation of Mann 0] and Wald 0. In particular, Wald introduced the notion of 
an algebra-valued tangent space (more generally, algebra-valued tensor fields) and gen- 
eralized the whole formalism of differential geometry (that is covariant derivatives. 
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metrics, forms and curvature) to the algebra-valued case. This geometric interpreta- 
tion was motivated by the preceding work of Cutler and Wald , where they found a 
new type of gauge invariance for a collection of massless spin-2 fields, more precisely, 
a consistent theory describing the interaction of a collection of massless spin-2 fields. 
The infinitesimal form of the new gauge transformations of the algebra-valued metric 
looks precisely like the infinitesimal diffeomorphisms but with algebra-valued infinites- 
imal vector fields. The consistency condition simply says that the commutator of two 
infinitesimal gauge transformations is an infinitesimal gauge transformation. The main 
conclusion of Cutler and Wald was that the consistency conditions can only be satisfied 
for associative commutative algebras. Therefore, all constructions "diagonalize" (every- 
thing commutes) and the interacting theory of a collection of massless spin-2 fields is 
simply a sum of the usual Einsten-Hilbert actions for each field (no cross- interactions). 
This restricts significantly the possible physical applications of this approach. Further, 
Wald noticed that the algebra-valued vector fields do not generate diffeomorphisms of 
a real manifold. Instead, they are diffeomorphisms on so-called 'algebra manifold', in 
which the coordinates themselves are algebra-valued. The supersymmetric extension of 



this approach was continued in [^ 



Similar results concerning consistent cross-interactions of a collection of massless 
spin-2 fields were obtained in |l^, where the multi-graviton theories were analyzed 
(from a different point of view). The authors of this work considered a collection 
of massless spin-2 tensor fields such that: i) the Lagrangian contains no more than 
two derivatives of the fields, ii) the interactions can be continuously switched on, and 
iii) in the limit of no interaction, the Lagrangian reduces to the sum of Pauli-Fierz 
Lagrangians for each field. The authors of this paper consider the free limit condition 
described above to be the crucial one. They mention that Cutler and Wald did not 
analyze the extra conditions that must be imposed on the gauge symmetries coming 
from the free limit condition. They indicate that none of the examples of Cutler 
and Wald models have the correct free limit, since some of the gravitons come with 
the wrong sign and, therefore, the energy of the theory is unbounded from below. 
The free limit requirement gives an extra constraint, which implies that the algebra 
is not only commutative and associate but also symmetric. Therefore, the algebra 
has a trivial structure: it is simply the direct sum of one-dimensional ideals. This 
eliminates all cross-interactions. The only consistent deformation of the sum of Pauli- 
Fierz Lagrangians is the sum of Einstein-Hilbert Lagrangians (with a cosmological 
term) for each field. 

In our previous papers [[l^, [ill we proposed to describe gravity by a multiplet 
of tensor fields with the corresponding gauge symmetry incorporated in the model. 
We would like to stress from the very beginning that our approach is different from 
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the schemes studied by previous authors. Our main ingredient is the matrix-valued 
two-tensor field, so that the components of these tensor fields do not commute with 
each other, in general. Our algebra is associative but noncommutative. The other 
difference (related to the first one) is the form of the gauge transformations. We start 
from the very beginning with a real manifold with real-valued coordinates and the 
usual diffeomorphism group, so that there is no problem of defining the finite gauge 
transformations. Our gauge group is simply the product of diffeomorphisms and an 
internal group, say U{N). That is why we do not have the inconsistences studied by H, 



3 and |l^ and we can allow our algebra to be non-commutative. In our approach there 
is a collection of tensor spin-2 fields, but only one of them is massless, all other fields are 
massive. How exactly this happens depends on the details of the symmetry breaking 
mechanism etc, but since we have only one (the usual real-valued) diffeomorphism 
group, only one field is massless. The parameter of our gauge transformations is a real 
valued vector field, not the algebra- valued vector fields of p, H, |T5[ needed to describe 
multiple massless spin-2 fields, which transform independently of each other. A similar 



approach to matrix gravity was proposed recently by Chamseddine [131, who considered 
an SL{2N, C) SL{2N, C) gauge theory with a spontaneous breakdown of symmetry 
giving one massless spin-2 field and (2iV^ — 1) massive spin-2 fields. 

Our idea was simply to repeat the Einstein's analysis of the causal structure of 
space-time. We showed that the basic notions of general relativity are based on the 
geometrical interpretation of the wave equation that describes propagation of fields 
without internal structure, (in particular, light), that could transmit information in the 
spacetime. At the microscopic distances this role of the electromagnetic field could be 
played by other fields with some internal structure. That is why instead of a scalar 
wave equation we studied a hyperbolic system of linear second-order partial differential 
equations. This cardinally changes the standard geometric interpretation of General 
Relativity. Exactly in the same way as a scalar equation defines Riemannian geometry, 
a system of wave equations generates a more general picture, that we call Matrix 
Geometry. 

Morevover, we developed some kind of "matrix geometry" . We defined the matrix- 
valued metric, affine connection, torsion and curvature. We have also found a particular 
form of a compatibility condition that enabled us to explicitly express the connection in 
terms of the derivatives of the metric. By using the matrix curvature we constructed an 
action functional that: i) contains no more than two derivatives of the fields in such a 
way that the equations of motion contain no more than two derivatives, ii) is invariant 
under both diffeomorphisms and the gauge transformations, iii) reduces to the sum of 
the Einstein-Hilbert action and the Yang-Mills action in the weak deformation limit. 
Notice that we do not obtain multiple Einstein-Hilbert actions, but just one. As we 
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mentioned above there are several choices to make when generahzing the real-valued 
geometry to the matrix- valued one. In particular, the definition of the "measure", 
raising and lowering indices, definition of higher-order covariant derivatives, torsion 
constraints, and some others. As a result the action constructed in these papers is not 
unique. 

In the present paper we propose an additional physical principle (Low Energy Limit) 
to construct an action functional that has all of the above properties but is also unique 
(or natural). We simply require that the classical action of gravity is the semiclassical 
limit of the effective action of matter fields interacting with the gravitational field. This 
is nothing but the ideology of "induced gravity". It is well known that the usual 
Einstein-Hilbert action (with a cosmological term) is just the first two terms of the 
asymptotic expansion of the effective action, or, simply, the first two coefficients of the 
asymptotic expansion of the trace of the heat kernel for a certain invariant second- 
order self-adjoint elliptic (in Euclidean formulation) partial differential operator with 
a scalar leading symbol given by the Riemannian metric. That is why, in the present 
paper we propose to define the action of matrix gravity in a similar way as a linear 
combination of the first two coefficients of the trace of the heat kernel for a more general 
differential operator, with a matrix-valued leading symbol given by the matrix-valued 
metric. This way we avoid the constructions of matrix geometry, such as curvatures 
etc, which are non-unique. Although the language of differential geometry is certainly 
helpful, we simply do not need it here. In this approach we get the invariants not from 
the curvatures hut from the spectral invariants of a differential operator with a matrix- 
valued leading symbol. This is the main idea and the main difference of the present 
paper from our previous papers |T0, in|. 



The outline of the present paper is as follows. In Sect. 2.1 we introduce the 
necessary differential-geometric framework. We define the vector space of complex N- 
vectors and N x N matrices. In Sect. 2.2 we introduce the fields as densities. This is 
done to avoid the definition of the measure, which we do not have since we do not have 
the metric. In Sect. 2.3 the matrix-valued metric and the matrix-valued Hamiltonian 
are introduced and certain non-degeneracy conditions are imposed. We consider two 
cases: elliptic and hyperbolic metric, depending on the eigenvalues of the Hamiltonian. 
In Sect. 2.4 it is shown that the above construction is closely related to a collection of 
Finsler geometries and the Finsler metrics are defined. In Sect. 2.5 we consider matrix- 
valued p-forms and define the matrix-valued star operator * (similar to the Hodge star 
operator). In order to define this operator in an invariant way we also introduce an 
additional ingredient: a matrix-valued measure. However, this measure is not related 
to the matrix- valued metric in the usual way (since it is not unique in the matrix case). 
Rather it is considered as an additional degree of freedom. 
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In Sect. 3 we define the invariant differential operators of the first and the second 
order acting on matrix- valued p-forms. We define the exterior derivative, the coderiva- 
tive as well as their covariant versions in Sect. 3.1 and the gauge curvature in Sect 3.2. 
By using the covariant derivative and the coderivative we define invariant second-order 
differential operators (Laplacians) in Sect. 3.3 and invariant first-order differential op- 
erators (Dirac operators) in Sect 3.4. In Sect. 3.5 we present the main operator studied 
in this paper in local coordinates. This is used later in the paper. 

In Sect. 4 we list several particular examples of decomposition of the general 
matrix valued metric and introduce the deformation parameter. We consider both 
commutative and noncommutative cases and compare them with the approaches of 
previous authors. 

In Sect. 5 we briefly summarize the heat kernel approach to compute the spectral 
asymptotics, the zeta-function regularization, and the asymptotic expansion of the 
effective action. The heat kernel coefficients introduced here play the crucial role in 
future development. In Sect. 5.1 we compute the first two coefficients of the asymptotic 
expansion of the trace of the heat kernel for the non-Laplace type differential operator 
of general type (that is, a second-order operator with a non-scalar leading symbol). 
The Theorem 1 is one of the main results of the present paper. In Sect. 5.2 we use the 
results of the Theorem 1 to calculate the coefficients Aq and Ai in the commutative 
(scalar) case and confirm the well known results. This serves as a 'consistency check' of 
the calculations of Sect. 5.1 as well as shows how to make use of the general formulas 
of Theorem 1. 

In Sect. 6 the results for the heat kernel asymptotics are used to construct the 
action of "matrix gravity" via the induced gravity approach. The form of the action 
in the commutative limit is provided. It contains a metric and a scalar field (coming 
from the additional degree of freedom associated with the measure). 

In Sect. 7 the results are briefiy discussed and a list of open interesting problems 
(singularities, high-energy behavior, dark energy, confinement etc) is presented. 



2. Vector Bundles 

2.1 Vector Spaces 

To be precise, let V^ be a iV-dimensional complex vector space with a Hermitian inner 
product ( , ). As usual, the dual vector space V* is naturally identified with V, ip ^-* i^, 
by using the inner product: 

(^,^)=^(^), (2.1) 
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where (fjip & V and ip & V*. If the elements (/? of the vector space V are represented 
by A^-dimensional contravariant vectors ip = (v^^^), {A = 1, . . . , A^), and the elements 
of the dual space V* by covariant vectors ip = {iPb), then the inner product is defined 
via the Hermitian metric, E = (Ej^^), 

{^, ^) = ij^E^ = i^^*Ej,j,^^ = i,BV^ , (2.2) 

where 

E^ = E, i> = i)^E, (2.3) 

or 

^s = ^''*Ej,^ = {E^Jj^r. (2.4) 

Here and everywhere below we use the Einstein summation convention and follow the 
standard notation to mark the indices of the complex conjugated fields by dots; the 
symbol * denotes the complex conjugation and the symbol f denotes the Hermitian 
conjugation. 

The vector space End {V) of endomorphisms of the vector space V is isomorphic 
to the vector space Mat (A^, C) of complex N x N matrices X = (X^b)- The group 
of automorphisms Aut {V) of the vector space V is isomorphic to the general linear 
group GL{N, C) of complex nondegenerate N x N matrices. The adjoint X of an 
endomorphism X is defined by 

X = E-^X^E, (2.5) 

so that for any ip,ip E V 

{^,X^) = {X^,^) . (2.6) 

The metric determines the subgroup of unitary endomorphisms (or isometries) G 
of the group of the automorphisms Aut {V) of the vector space V preserving the inner 
product, that is 

G = {U e GL{N, C) I f/f/ = 1} , (2.7) 

where I is the identity endomorphism. The dimension of the group G is dimG = A^^. 
In the case E = I the group G is nothing but U{N). 

2.2 Vector Bundles of Densities 

Now, let M be a smooth compact orientable n-dimensional manifold without boundary 
and V [w] be a smooth vector bundle of densities of weight w over the manifold M 
with the fiber V. Here and everywhere below we indicate explicitly the weight in the 
notation of the vector bundles. That is, the sections ip of the vector bundle V[w] 
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are vector- valued functions '^{x) = (ip'^^x)) that transform under the diffeomorphisms 
x'^ = x'^{x) according to 

^'{x') = J~'"{xMx), (2.8) 



where 

J{x) = det ^ ^ 



dx'^ 



(2.9) 



Hereafter we label the local coordinates on the manifold M by Greek indices which run 
over 0, 1, . . . , n — 1. The small Latin indices will be used to label the space coordinates, 
i.e. the local coordinates on hypersurfaces transversal to the time coordinate; they will 
run over 1, 2, . . . , n — 1. 

Further, we denote by V* [w\ the dual bundle of densities of weight w with the fiber 
V* and by End (V) [w\ the bundle of endomorphisms-valued densities of weight w of 
the bundle V [w] with elements being matrix- valued functions U = {U^b{x)). 

We will also consider the bundle of vector- valued antisymmetric densities of weight 
w of type (0,p) (p-forms), Ap[w] = ((A^T*M) (g) V) [w], and the bundle of anti-symmetric 
densities of weight w of type (p, 0), A^lw] = ((A^TM) ® V) [w]. Finally, we also in- 
troduce the following notation Ep[w] = ((A*'T*M) End(V)) [w] for the bundle of 
endomorphism-valued densities of weight w of type (0,p) and Ep[w] = {{A^TM)(E) 
End (V))[w] for the bundle of endomorphism-valued densities of weight w of type {p, 0). 

The group G is also promoted to a smooth vector bundle ^[0] (gauge group), a 
subbundle of the bundle End(V)[0]. Under the action of the gauge group ^[0] the 
sections if of the bundles V[w] and V*[t(7] transform as 

ip'{x) = U{x)ip{x) , ^'(x) = >^{x)U-\x) , (2.10) 

The metric E is assumed to be a section of the bundle (V* ® V*) [a] with some 
weight a. It is invariant under the gauge transformations 

E'{x) = {U\x))-^E{x)U-\x) = E{x) , (2.11) 

which guarantees the invariance of the fiber inner product. 

However, to get a diffeomorphism-invariant L^ inner product 

{^,V) = Jdx{^{x),<f{x)) , (2.12) 

M 

and the L^ norm 

II<^IP = i'P,^) = dx {ip{x),ip{x)) , (2.13) 

M 



on the vector bundle V[w] with the standard Lebesgue measure dx = dx^ A ■ ■ ■ A rfx", 
the metric must be a density of weight (1 — 2w), more precisely, a section of the bundle 
(V* (S> V*)[l — 2w]. Then the completion of C°°(V[w]) in this norm defines the Hilbert 
space L'^{V[w]). Notice that this means, that on the bundle V[|] of densities of weight 
w = ^ the metric E{x) is invariant under the diffeomorphisms as well, and, therefore, 
without loss of generality can be assumed to be constant. 

2.3 Noncommutative Metric 

Our goal is to construct covariant self-adjoint second order differential operators acting 
on a smooth sections of the bundles Ap[|] and A^[i], that are covariant under both 
diffeomorphisms , 

LV'(x') = J-^'^{x)L^{x) , (2.14) 

and the gauge transformations 

L'^\x) = U{x)Lip{x) . (2.15) 

To do this we need to define the following objects. First of all, we need an isomor- 
phism between the bundles Ap[|] and Ap[|]. This is usually achieved by the Hodge star 
operator, which is defined with the help of a Riemannian metric. This is exactly the 
place where we want to generalize the standard theory since we do not want to introduce 
a Riemannian metric; instead of a Riemannian metric (which is an isomorphism be- 
tween the tangent TM and the cotangent T*M bundles) we introduce an isomorphism 
between the bundles Ai[w] = {T*M O V)H and K^[w] = {TM V)[w], i.e. 

a:AiH^A^H. (2.16) 

Such an isomorphism is determined by a section of the vector bundle [TM TM 
End(V))[0] defined by the matrix-valued symmetric tensor a'^'^ = [a^'^^Bix)) that 
satisfies the following conditions: 

1. The matrix a''" is self-adjoint (recall that aF" = E~^{a^"'yE), 

W = a^" (2.17) 

and symmetric, 

a"^ = a'^^. (2.18) 

2. It transforms under the diffeomorphisms as 

a'^'^(x') = ^^a"^(x), (2.19) 
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and under the gauge transformations as 

a'''^(a;) = U{x)a>"'{x)U-\x) . (2.20) 

3. Consider the matrix 

H{x,0=a^''{x)U^'^ (2-21) 

with ^ being a cotangent vector ,^ at a spacetime point x. Since this matrix is 
self-adjoint, all its eigenvalues /i(j)(x,^), (i = l,...,s), must be real. We will 
assume that they have constant multiplicities di. The eigenvalues are invariant 
under the gauge transformations ( p.lOD and transform under the diffeomorphisms 
as 

K{x',0 = Hx,e), (2.22) 

where 

We will also require that this matrix satisfies one of the following conditions. 
Elliptic case: The matrix H{x,^) is positive definite, i.e. all its eigenvalues are 
positive, /i(j)(x,^) > 0, for any x and any ^ 7^ 0. Hyperbolic case: There is a 
one-form r = r(x) (specifying the time direction) such that at each point x for 
each eigenvalue /i(j)(a;,r(x)) < 0, and for any cotangent vector ( ^ not parallel 
to r each equation /i(j)(x, C + Ar(x)) = 0, has exactly two real distinct roots, 
A = A±(x,C). More precisely, a second-order partial differential operator L with 
the principal symbol H{x, ^) is hyperbolic if the roots Xj{x, () of the characteristic 
equation 

det H{x, C + Ar(x)) = (2.24) 



are real and strictly hyperbolic if they are real and distinct p5|, p9|). 



The infinitesimal forms of the diffeomorphisms and the gauge transformations are 
obtained as follows. Let x' = ft{x) be a one-parameter subgroup of diffeomorphisms 
such that 

ft{x)\t=o = x, and ^{x) = ^ft{x) . (2.25) 

at t=o 

Then 

S^a"" = ^a'^^{x') = Cdaa''^ - (<9„e'')«°'' - (9^0^"'" , (2.26) 

Similarly, let Ut be a one-parameter subgroup of the gauge group such that 



d 

Ut\t=o = I, and u = —Ut 

dt 
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(2.27) 
t=o 



Then 

at 



= [uj, an . (2.28) 



Notice that these infinitesimal transformations are different from the gauge trans- 
formations studied in P, H, |l5l since there is only one vector field ^. 



2.4 Finsler Geometry 



The above construction is closely related to Finsler geometry |28]. First of all, we note 
that the eigenvalues h(^i){x,^) are homogeneous functions of ^ of degree 2, i.e. 

%)(x,A0 = A2%)(x,0. (2.29) 

Next, we define the Finsler metrics 

<.)(^,0 = ^^|^/^W(^>0- (2.30) 

All these metrics are non-degenerate. In the elliptic case all metrics g'f^^ are positive 
definite; in the hyperbolic case they have the signature ( — !-■■■+). In the case when 
a Finsler metric does not depend on ^ it is simply a Riemannian metric. 
The Finsler metrics are homogeneous functions of ^ of degree 

9^{x,X0 = 9^ix,0, (2.31) 

so that they depend only on the direction of the covector ^ but not on its magnitude. 
This leads to a number of identities, in particular, 

hi^){x,0=9'^:^{x,0U^ (2.32) 

and 

^^h^,^{x,0=9^;^{x,0^u- (2.33) 

Next we define the inverse (covariant) Finsler metrics by 

9i^),.ix,x)g'(^ix,0=S';, (2.34) 

where x'^ is the tangent vector defined by 

i;'^ = ^[:j(x,Oe., (2.35) 

so that 

^f^ = 9ii^^iuix,x)x'' . (2.36) 
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Finally, this enables one to define the Finsler intervals 

dsfi) = g{i),iu{x, x)dx^dx'' . (2.37) 

The existense of the Finsler metrics allows one to define various connections, curvatures 
etc (for details see p8[| ). 

As we see, the propagation of gauge fields induces Finsler geometry. 

2.5 Star Operators 

Since the map a (|2.16|) is an isomorphism, the inverse map 

6 = a-^ : A^H ^ AiH. (2.38) 

is well defined. In other words, for any ip^ = (ip^^^) there is a unique ipi, = (v^^) 
satisfying the equation a^"ipy = ip^, and, therefore, there is a unique solution of the 
equations 

a^''K^ = 6^, b^,a''^ = Si^. (2.39) 

Notice that the matrix b^^ has the property 

^M- = b,^ , (2.40) 

but is neither symmetric 6^,^ ^ h^^ nor self-adjoint h^^ ^ h^y. 

The isomorphism a naturally defines the maps between the bundles Ap[w] and 

A : ApH -^ KP[w] , B : A^H -^ Kp[w] , (2.41) 



as follows 
where 
and 
where 



{A^Y^■■■^'^ = A'^i-'^^^i-^''(/?^^...^^ , (2.42) 

^Mi-Mp.i-.p = Alt ^,...^^ Alt ,,...,y' ■ ■ ■ a'''"- (2.43) 

{Bip),,...^^ = B^,...,^,,...,y--''-, (2.44) 



Here Alt^^.-.^p denotes the complete antisymmetrization over the indices /xi, . . . ,/ip. 

We will assume that these maps are isomorphisms as well. Strictly speaking, one 
has to prove this. This is certainly true for the weakly deformed maps (maps close to 
the identity; more on this later). Then the inverse operator 

A-^ : A^[w] ^ ApH , (2.46) 
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is defined by 

(A-V),....,, = A-l.^^,^...,^^^--^^ , (2.47) 

where A"^ is determined by the equation 

A-} „ , , ^^i-^p"i-a^ = <5"^ ■ ■ ■ r^, . (2.48) 

Notice that because of the noncommutativity, the inverse operator A~^ is not equal to 
the operator B, so that A~^B ^ Id. 

This is used further to define the natural fiber inner product on the space of p-forms 

Ap via 

(^, V) = ^^^,...^, A^^-^-'^^-'^-V^.,..., . (2.49) 

We will also need a smooth self-adjoint section p (given locally by the matrix- valued 
function p = (p'^six))) of the bundle End (V)[|] of endomorphism- valued densities of 
weight J. That is p satisfies the equation 

p = p, (2.50) 

(recall that p = E~^p'^E), and transforms under diffeomorphisms as 

p'{x') = J-'/\x)p{x) , (2.51) 

and under the action of the gauge group ^[0] as 

p{x) = U{x)p{x)U-\x) . (2.52) 

Clearly, the matrices 

a^^ = g^^l, p = g'l'\. (2.53) 

where g^^ is a (pseudo)-Riemannian metric and 

c/= |detc/'^n"\ (2-54) 

satisfy all the above conditions. We will refer to this particular case as the commutative 
limit. 

Of course, (on orientable manifolds) we always have the standard volume form e, 
which is a section of the bundle -£'„[—!] given by the completely antisymmetric Levi- 
Civita symbol Sm-'-n^- The contravariant Levi-Civita symbol e with components 

e^i-Mn ^ ^^^^^^^^^ ^ (2.55) 
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with a = +1 in the elhptic case and o" = — 1 in the hyperbohc case, is a section of the 
bundle E"[l]. 

These densities are used to define the standard isomorphisms between the densities 
bundles 

e : AP[w] ^ An-p[w - 1], £: ApH ^A"-P[w + 1] (2.56) 

by 

(^V^),....,„_, = ^e,,...,„_^.,....y--^ , {e^r--^-^ = l^A.--M"-.-i-.v^,^...,^ . (2.57) 

By using the well known identity 

e,,...,„_,.,....,e''^-^-^'^-'^ = a{n-p)\p\6l\ ■ ■ -5^;] (2.58) 

we get 



ee 



e£ = a(-l)p(«-p)ld . (2.59) 



By combining e and e with the endomorphism p we can get the invariant form 
ep"^, which is a section of the bundle -E„[0], and the form ep~'^, a section of the bundle 
i5"[0]. Notice, however, that, in general, the contravariant form ep^'^ is not equal to 
that obtained by raising indices of the covariant form ep^, i.e. £p~^ 7^ Aep^ or 

If we require this to be the case then the matrix p should be defined by 

p = V-'^\ (2.61) 

where 

V = <y-e^^...^^e,^...,^a^-^^ ■ ■ ■ a^-'- . (2.62) 

Clearly, 77 is a density of weight (—2). Since a^^ is self-adjoint, we also find that 77 and, 
hence, p is self-adjoint. The problem is that in general 77 is not positive definite. Notice 
that in the commutative limit 

r7 = (Tdet(7^"=|det(7^,|-\ (2.63) 

which is strictly positive. 

Therefore, we can finally define two different star operators 

*,S: ApH -^ A„_pM (2.64) 
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by 

* = epAp , S = p-^A-^p-^e (2.65) 

that is 

(*^)m-M„-p = ^^«-Mn-p^i-pP^'^"''^^""''"VV^ar-ap , (2.66) 

The star operators are self-adjoint in the sense 

(<^, *V^) = (*^, ^) , (<^, *^) = (*<^, ^) , (2.68) 

and satisfy the relation: for any p form 

*i = S* = a(_i)p("-p) Id . (2.69) 

2.6 Hilbert Spaces 

Now, let us consider the bundles of densities of weight |, V[|], and, more generally, 
Ap[i]. If dx = dx^ A ■ ■ ■ A dx" is the standard Lebesgue measure in a local chart on M, 
then we define the diffeomorphism-invariant L^-inner product 



(^,V9) = y rfx(^(x),^(a;)) , (2.70) 

M 

and the L^ norm 

||<^||2=(v9,v9) = y"rfx(v9(x),v9(x)) , (2.71) 

M 

The completion of C°°(Ap[i]) in this norm defines the Hilbert space L^(Ap[|]). 

3. Differential Operators 

3.1 Noncommutative Exterior Calculus 

Next, we define the invariant differential operators on smooth sections of the bundles 
Ap[0] and A*'[l]. The exterior derivative (the gradient) 

d:C°°(Ap[0])^C°°(Ap+i[0]) (3.1) 

is defined by 

(c^V^)mi-mp+i = (p+1)%iV'm2-mp]' if p = 0,...,n-l, (3.2) 
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dip = a p = n , (3.3) 

where the square brackets denote the complete antisynimetrization. The coderivative 
(the divergence) 

d : C°"(AP[1]) ^ C^{A^-^[1]) (3.4) 

is defined by 

d = ai-ir^+^ede . (3.5) 



By using ( |2.58|) one can easily find 

{dipY'->'^-' = S^^'^^i-'^^-i if p = 1, . . . , n , (3.6) 

difi = if p = . (3.7) 

One can also show that these definitions are covariant and satisfy the standard relations 

d^ = (P = 0. (3.8) 

The endomorphisni p is a section of the bundle End(V)[|]. Therefore, if y) is a 
section of the bundle Ap[i], the quantity p~V is a section of the bundle Ap[0]. Hence, 
the derivative d{p~^{p) is well defined as a smooth section of the vector bundle Ap+i[0]. 
By scaling back with the factor p we get an invariant differential operator 

pdp-': C°°(A,[i])^C°°(A,+i[i]). (3.9) 

Similarly, we can define the invariant operator of codifferentiation on densities of weight 
1 

p-'dp : C^ (A^ [i]) -. C^ (A^-i [i]) . (3.10) 

Now, let ;B be a smooth anti-self-adjoint section of the vector bundle £^i[0], defined 
by the matrix- valued covector B^ = {Bi^i^b{x)) i.e. 

B^ = -B^, (3.11) 

that transforms under diffeomorphisms as 

^;(^') = ^^.(^), (3.12) 

and under the gauge transformations as 

B'^ix) = U{x)B^{x)U-\x) - {d^U{x))U-\x) . (3.13) 

Such a section naturally defines the maps: 

i3:ApH^Ap+iH (3.14) 
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by 

(%)/.i-mp+i = (p + l)^[m¥'M2-Mp+i] (3-15) 

and 

B : KP[w] -^ M-\w] (3.16) 

by 

Notice that 

B = (T{-l)''P+^iBe (3.18) 

similar to (|375|) . 

This enables us to define the covariant exterior derivative 

P:C°°(A,[i])->C-(A,+i[i]). (3.19) 

by 

V = p{d + B)p-^ (3.20) 

and the covariant coderivative 

V : C°° {Af [I]) ^ C°° {AP'^ [I]) , (3.21) 

by 

V = a{-l)''P+^iVe = p-\d + B)p. (3.22) 

These operators transform covariantly under both the diffeomorphisms and the gauge 
t r ansf ormat ions . 

3.2 Noncommutative Gauge Curvature 

One can easily show that the square of the operators V and f> 

p2:C-(A,[i])^C~(A,+2[i]) (3.23) 

f)2 . C^^f^P+2 [1]) ^ ^oo^^p ^ij) (3_24) 

are zero-order differential operators. In particular, in the case p = they define the 
gauge curvature JF, which is a section of the bundle -E2[0], by 

(PV)m- = pJ^f.up-^f , ^V = P'^^F^uP'p''^ , (3.25) 

where 

J^ = dB + [B, B] , (3.26) 
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and the brackets [ , ] denote the Lie bracket of two matrix-valued 1-forms, i.e. 

[A, B]^, = A^B, - B,A^ . (3.27) 

The gauge curvature is anti-self-adjoint 

^fiv = —^fiu (3.28) 

and transforms covariantly under diffeomorphisms 

and under the gauge transformations 

^lui^) = U{x)T^u{x)U-\x) . (3.30) 

3.3 Noncommutative Laplacians 

Finally, by using the objects introduced above we can define second-order differential 
operators that are covariant under both diffeomorphisms, and the gauge transforma- 
tions. In order to do that we need first-order differential operators (divergences) 

V: C-(A,[i])^C°°(A,_i[i]), (3.31) 

First of all, by using the L^ inner product on the bundle Ap[|] (recall that the 
metric E is constant in this case) we define the adjoint operator V by 

{if^Vi,) = {V^.i,) . (3.32) 

This gives 

V = -A-^VA = -a{-l)''P+^A~^iVeA = -A'^p-^d + B)pA , (3.33) 

which in local coordinates reads 

(V^),,...,, = -{A-'),,...,^,,....^p-\d, + S,)pA— ^— -+Vai...a,+i , (3.34) 

The problem with this definition is that usually it is difficult to find the matrix 

[A )lJ.l---IJ.pUl---Up- 

Then we define the second order operators 

VV ,VV,A: C^ (A, [i]) ^ C^ (A, [i]) , (3.35) 
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where the "noncommutative Laplacian" is 

A = -VV - VV 

= A-^VAV + VA-^VA 

= A-^p-\d + B)pAp{d + B)p-^ 

+p{d + l3)p-^A-^p-\d + B)pA. (3.36) 

In local coordinates this reads 

(A^)^,..^^ = |(p + l)A-\.^,.,....,p-'(9. + S.)p 
xA'"''-''^"'''-"^p{d^ + B^)p-^ 

xp~\d, + B,)pA^^--^---^--^A^^„„^^ (3.37) 

We could have also defined the coderivatives by 

Vi = -*V*, f)2 = -BVA , "Dg = -*V* , f)i = -*V*. (3.38) 

These operators have the advantage that "Di is polynomial in the matrix a^^ and 'D2 
is polynomial in the matrices a^^ and h^^. However, the second order operators "DjI), 
'D'Dj and Aj = —'DjV — 'D'Dj, {j = 1,2,3,4), are not self-adjoint, in general. In the 
commutative limit all these definitions coincide with the standard de Rham Laplacian. 
In the special case p = the "noncommutative Laplacian" A reads 

A = p-\d + B)pAp{d + B)p-^ , (3.39) 

which in local coordinates has the form 

A = p-\d^ + B,)pa^^p{d, + B,)p-^ . (3.40) 

3.4 Noncommutative Dirac Operator 

It is worth mentioning another approach. Suppose we are given a self-adjoint section 
r of the bundle -Ei[0], given locally by matrix valued vector V^ = (r^^B(a:)) satisfying 

^^ = 1^, (3.41) 

where as usual T^ = E'^lV^yE. It transforms under the diffeomorphisms as 

r"^(^') = ^r'^(x) (3.42) 
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and under the gauge transformation as 

r'^(a;) = U{x)r^'{x)U-\x) . (3.43) 

Then the matrix 

^^^u ^ _ (p^r^ + r^r^) (3,44) 

is self-adjoint a'^'^ = a^^'^ and symmetric a'^'^ = d'^^. Moreover, suppose this matrix 
satisfies all the conditions for the matrix a listed above in Sect. 2.3. Most importantly, 
the matrix 

H{x, i) = a>^''{x)i,i, = [V{x, i)]\ (3.45) 

where 

V{x,0='^^{x)i,, (3.46) 

is nondegenerate and satisfies the ellipticity or hyperbolicity conditions formulated in 
Sect. 2.3. (The advantage of this approach is that in the elliptic case the matrix H{x, ^) 
is automatically positive definite.) 

Such vector naturally defines the map 



by 

and the map 

by 



r : C^{KP[w]) ^ C°°(AP+^H) (3.47) 

^Y^Y'-^''+' = {p+ l)r[^V''''"^''+'' (3.48) 

f : C"-(ApM) ^ C"-(Ap_iH) (3.49) 



(rv^)M....M.-i = r^v^^m...Mp • (3.50) 

Therefore, we can define first-order invariant differential operator ( "noncommuta- 
tive Dirac operator") 

D:C-(A,[i])^C°°(A,[i]) (3.51) 

by 

D = VV = Vp{d + B)p~^ , (3.52) 

which in local coordinates reads 

{D^),,...,, = {p + l)r^p(5[^ + %)P^ Vmi...m.] • (3.53) 

The adjoint of this operator with respect to the L^ inner product is 

D = -A-^VTA = -A-^p-\d + B)prA , (3.54) 
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which in local coordinates becomes 

(^V^),,...,, = -ip + lX...,,.,....,p-'(5. + S.)pr['^A'^-->--^v^.,...„^ . (3.55) 

These can be used then to define the second order self-adjoint differential operators 
( "noncommutative Laplacians" ) 

Ai = -DD 

= TVA-^VTA 

= rp{d + B)p-^A-^p-\d + B)prA, (3.56) 

A2 = -DD 

= A-^VTATV 

= A-^p-\d + B)pTAtp{d + B)p-^ . (3.57) 

In the case p = these operators have the form 

Ai = VVVT 

= T^p{d^ + B^)p'\d, + B,)pV- , (3.58) 

A2 = VTTV 

= p-\d, + B,)pV-V>^p{d, + B^)p-^ . (3.59) 

In the elliptic case the above constructions can be used to develop noncommutative 
generalization of the standard Hodge-de Rham theory, in particular, noncommutative 
versions of the index theorems, the cohomology groups, the heat kernel etc. This is a 
very interesting topic that requires further study. 

3.5 Differential Operators for p = 

In present paper we will restrict ourselves to the case p = 0, more precisely, to the 
second order self-adjoint differential operators acting on smooth sections of the bundle 
V[|] that can be presented in the form 

L = -p-\d^ + B^)pa^'p{d, + B,)p-^ + Q , (3.60) 

by an appropriate choice of a^^ , p, B^ and Q G C°°( End (V)[0]). It is this operator 
that we will study in detail below. 

This operator can be also written in the form 

L = XX"^z. + Q, (3.61) 
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where 

X^ = d^ + C^, (3.62) 

X^ = -d^ + C^ , (3.63) 

and 

C^ = -p,^p-^+pB^p-\ (3.64) 

Cf, = -p'^p^L - p'^B^p , (3.65) 

comma denotes partial derivatives, e.g. p_^ = d^p, 
In more details, it can be presented in the form 

L = -a^'d^d, + Yf^d^ + Z , (3.66) 

where 

Y^" = -a"",^ + Cyd"" - a^T^ (3.67) 

Z = Q- (a'^T,),^ + C^a^T, (3.68) 

The leading symbol of the operator L is given by the matrix -ff (x, ^) = ci^'^{x)C,iiC,u, 
with ^ a cotangent vector. In the Sect 2.3 we formulated precise conditions on this 
matrix for the operator L to be elliptic or hyperbolic. The system of hyperbolic partial 
differential equations describes the propagation of a collection of waves and generates 



the causal structure on the spacetime manifold |T^, |TT|. In the following we restrict 
ourselves for simplicity to the elliptic case, in which the leading symbol H{x,C,) is 
positive definite. 

4. Examples 

We see that in the matrix case the operator L does not define a unique Riemannian 
metric g'^'^. Rather there is a matrix- valued symmetric 2-tensor field a'^'^ . Although the 
matrix g^'^ plays the role of the Riemannian metric in the commutative limit, it loses 
this role in fully noncommutative strongly deformed theory. 

Abelian Case. This is the case studied by Cutler and Wald ||, § and Boulanger et 



al |T5[. In this case the matrix a^'^ has the following form 

a'^'^ = <)%) , P = ^;iX), (4.1) 

where g(^a) = ( det g'^^))~^, and 11(a) are the projections on the a-th component, that 

\a)^ 



is diagonal matrices such that {I[(^a))^B = ^B^fa)^B (^° summation !). In this case 
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all the structures commute and there is really nothing new in the matrix gravity 
theory, which is just the sum of the Einsten-Hilbert actions. What we consider 
in the present paper is radically different since it involves fully noncommutative 
theory. Some of the examples are considered below. 

Commutative Limit. The simplest case is when there is a decomposition 

a^^ = g^'n + x/i^'' , p = g"!^ exp(x0) , (4.2) 

where I is the unity matrix, g^^ is a contravariant symmetric 2-tensor, h}^^ is a 
matrix-valued tensor field, is a matrix- valued scalar field and x is a deformation 
parameter. 

Reducible Case. A more general decomposition is given by 

a^"^ = g^Tl + xh^" , p = g^'^n-'"'"^ exp(x0)fi-"/s , (4.3) 

where f2 is a self-adjoint positive definite (non-constant) section of the endomor- 
phism bundle End (V)[0], i.e. H = fi, f] > 0. 

Commutative N = 2 Case. The simplest nontrivial model is when the matrix a^'^ is 
a real symmetric 2x2 matrix with the gauge group 0(2) and the fiber metric 
E = 1. We have then the following decomposition 

af"' = gf"'! + >iTh^"' , p = g^l^ exp(xr0) , (4.4) 

where h}^^ is a tensor field, is a scalar and 

The condition of positivity of the matrix H[x^ ^) = a^^^^^v reads 

>c\h^''UA < \9'''UA (4.6) 

for any .^ 7^ 0. 

Noncommutative N = 2 Case. The simplest example of a fully deformed noncom- 
mutative theory is the case of A^ = 2 with the gauge group U{2). Assuming that 
the fiber metric is given by 

E = I, (4.7) 
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we have the decomposition in terms of Pauh matrices 

a^- = g^^-l + xr./if , p = g^'^ exp(xr„</)J (4.8) 

where h'^^ are tensor fields, 0^ are scalars and r^ are standard Pauh matrices 

-(°i)' -C"o)' -(o-O- <-' 

The condition of positivity of the matrix H{x,C,) = a'^'^^^iC.u reads 

{g^'^g^^ - x^/if/if ) UuUp > (4.10) 

for any ^ 7^ 0. 

Large A^ Case. It is also interesting to consider the case of the U{N) gauge group 
in the limit A^ ^ 00. Since the matrix- valued metrics do not commute, it could 
serve as a toy model for quantum gravity. The decomposition has the form 

a^"" = g^"'I + x{N)h'^''Ta, (4.11) 

where T^ are the generators of SU{N) and one can adjust the dependence of the 
coupling X on iV to simplify the limit N ^ 00. 

5. Spectral Asymptotics 

As we mentioned above we restrict ourselves to elliptic operators. All the geometric 
quantities we need for the gravity theory can be extracted from the spectrum of the 
elliptic operator L (p.60 ) (on a compact manifold M). The spectrum of the operator 



L is defined as usual by 

Lipk = Xkffk, llv^fclUa = 1, (5.1) 

where k = 1,2, . . .. It is well known that a self-adjoint elliptic partial differential 
operator with positive definite leading symbol on a compact manifold without boundary 



has a discrete real spectrum bounded from below p6| . Since the operator L transforms 
covariantly under the diffeomorphisms as well as under the gauge transformations (p.lO|) 
the spectrum is invariant under these transformations. 

To study the spectral asymptotics, i.e. the behavior of Afc as A; — * cxd, one introduces 
the spectral functions, in particular, the heat trace 



00 



Q{t) = y e-''\ (5.2) 



fc=i 
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and the zeta function 

oo 

C(.) = /i^^5^(A, + m2)-% (5.3) 

fc=i 

where each eigenvalue is taken with its multiphcity and /x is a renormahzation parameter 
introduced to preserve dimensions. Here t is a real parameter, m^ is a sufficiently 
large positive mass parameter (so that the operator {L + m^) is positive) and s is a 
complex parameter. It is well known that for a self-adjoint elliptic second-order partial 
differential operator with positive definite leading symbol the series ( ^.2| ) converges 
and defines a smooth function for t > and the series ( ^.3|) converges for Re s > n/2 
and defines a meromorphic function with simple isolated poles on the real line. These 
spectral functions are related by the Laplace-Mellin transform 

oo 



as) = ^J dtf-'e''^'e{t). (5.4) 



Moreover, it is also known that ({s) is analytic at s = 0, which enables one to define 
the regularized determinant of the operator {F + rn^) (the one-loop effective action in 
quantum theory) [Q 

1, _ fL + m'^\ 1 



2 ° \ K J 2 



r(i) = - log Det — = --d.Cis) . (5.5) 



s=0 



It turns out that the study of the spectral asymptotics is equivalent to the study of 
the asymptotic expansion of the heat trace 0(t) as t — *^ and that of the zeta function 
and the effective action as m^ — ;► oo. It is well known that there is an asymptotic 



expansion of the heat trace invariant as t ^ 0^ p6 



e(t)~(4vr)-"/2^t'=-iA,, (5.6) 

fc=0 

and an asymptotic expansion of the zeta-function as m — > oo 

C(.) ~ (47r)-/2 Yl rn " f'"'^'"''-"'^^ ' (5-7) 

fc=0 ^^' 



where Ak are so-called heat invariants 26, 20, 16, 19, 18 . The coefficients Ak are 



spectral invariants determined by the integrals over the manifold of some local invariants 
[^ |l^, |20| (for a review, see [|^, |2l]]) constructed from the coefficients of the operator 
L and their derivatives so that they are polynomial in the derivatives of the coefficients 
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of the operator L. The heat invariants Ak determine further the large mass asymptotic 



expansion of the effective action as m — > oo [16, 19, 20, 21 



r(i) ~ Com"y4o + Cim""^y4i + ^ CkivJ'"'^^ A^ 



k=2 



where 



Ck = --(47r) 



-n/2 



d_ 
ds 



T{s) 



m 



/i 



-2s 



(5. 



(5.9) 



s=0 



One of the most powerful methods to study the spectral asymptotics is the heat 
kernel. The heat kernel, f/(t|a;, x'), is defined as the kernel of the heat semigroup 
exp(— tL) operator for t > 0, i.e. the fundamental solution of the heat equation 



with the initial condition 



{dt + L)U{t\x,x') =0 



f/(0+|a;,x') =5{x-x') 



(5.10) 
(5.11) 



5{x — x') being the Dirac distribution. 

For t > the heat kernel U{t\x, x') is a smooth function near the diagonal of M x M 
and has a well defined diagonal value 



f/(t|a:,a;)=^e- 



t\i 



(pk[x) (S)(pk[x) 



(5.12) 



fc=i 



which has the asymptotic expansion as t — ;> p6 



U{t\x,x) ~ (47r)-'^/2^t^-tafc(x) 

fc=0 



(5.13) 



where ak{x) are some covariant densities. 

Moreover, the heat semigroup exp(— tL) is a trace-class operator with a well defined 
L^ trace 

Tri2exp(-tL) = dx tTvU{t\x,x) = Q{t) , (5.14) 

M 

where tr y denotes the fiber trace. We have defined the heat kernel in such a way 
that it transforms as a density of weight | at both points x and x'. More precisely, 
it is a section of the exterior tensor product bundle V[|] Kl V*[|]. Therefore, the heat 
kernel diagonal transforms as a density of weight 1, i.e. it is a section of the bundle 
End(V)[l], and the trace Tr ^2 exp{—tL) is invariant under diffeomorphisms. 
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Therefore, the heat invariants Ak can be constructed by computing the t ^ 
asymptotics of the solution of the heat equation and integrating the coefficients ak{x) 
over the manifold, i.e. 

Ak = dx tivo-kix) . (5.15) 

M 

A second-order differential operator is called Laplace type if it has a scalar leading 
symbol. Most of the calculations in quantum field theory and spectral geometry are 
restricted to the Laplace type operators for which nice theory of heat kernel asymptotics 



is available [^ |T6|, [1^, |1^, ^, |21|. However, the operators condidered in the present 
paper have a matrix valued principal symbol H{x, ^) and are, therefore, not of Laplace 
type. The study of heat kernel asymptotics for non-Laplace type operators is quite new 
and the methodology is still underdeveloped. As a result even the first heat invariants 



{Aq, Ai and A2) are not known in general. For some partial results see |23, |22|, ^ 

5.1 Calculation of the Invariants Aq and Ai 

For so called natural non-Laplace type differential operators, which are constructed 
from a Riemannian metric and canonical connections on spin-tensor bundles the coef- 
ficients Ao and Ai were computed in |2^. Following this paper we will use a formal 



method that is sufficient for our purposes of computing the asymptotics of the heat 
trace of the second-order elliptic self-adjoint operator L ( p.60| ). 



First, we present the heat kernel diagonal in the form 

U{t\x, x) = I j^e-^^^ exp(-tL)e^«- , (5.16) 

where ^x = ^fj,x^, which can be transformed to 

Uit\x, x) = j -^ exp [-t iH + K + L)]-I, (5.17) 

R" 

where H is the leading symbol of the operator L 

H = a'''^^^,, (5.18) 

and K is a first-order self-adjoint operator defined by 

K = li^ [X.a^" - a'^'X,) , (5.19) 



where the operator X^ is defined by ( |3.62|) . Here the operators in the exponent act on 



the unity matrix I from the left. By changing the integration variable ^ ^ t ^^^^ we 
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obtain 



U{t\x,x) = (47rt)-"/2 



/ ^ ^^P (-^ -ViK-tiyi, (5.20) 



and the problem becomes now to evaluate the first three terms of the asymptotic 
expansion of this integral as t ^ 0. By using the Volterra series 



exp{A + B) = e^ + ^ dTk dru-i ■■■ I dn 







(5.21) 



we get 



exp (-H - VtK - tL\ = e-^ 
1 



_ 1 T2 1 





2\ 



+0{t') . 



(5.22) 



Now, since K is linear in ^ the term proportional to t^^^ vanishes after integration over 
^. Thus, we obtain the first two coefficients of the heat kernel diagonal 



U{t\x, x) ~ (47rt)-"/2 [ao + tai + 0{t^)] 



in the form 



oo 



d^,^-H 



TT 



n/2 



(5.23) 



(5.24) 



1 T2 



a,= f-^ /"rfr2 /"rfrie-(^--)^ire-(---)^ire--^^ 










(5.25) 
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These quantities are matrix densities that define finally the heat invariants 



An 



dx 



TX 



n/2 



trye 



-H 



(5.26) 



M 



i T2 



M 



i 



(5.27) 



These quantities are invariant under both the diffeomorphisms and the gauge transfor- 
mations. The coefficient Qq is constructed from the matrix a but not its derivatives, 
whereas the coefficient ai is constructed from the matrix a and its first and second 
derivatives as well as from the first derivatives of the field B and the matrix p and its 
first and second derivatives; obviously, it is linear in the matrix Q. Morevover, it does 
not depend on the derivatives of Q and is polynomial in the derivatives of a^*^, p and 
i3^, more precisely, linear in second derivatives of a and p and the first derivatives of 
B and quadratic in first derivatives of a and p. Further, since the operator L is self- 
adjoint, the heat kernel diagonal is self-adjoint matrix density and the heat trace is a 
real invariant. Therefore, the coefficients ao and ai are self-adjoint, and the invariants 
Aq and Ai are real. 

By integrating by parts and using the cyclic property of the trace one can simplify 
the expression for the invariant Ai as follows 

1 



A, 



dx 



d^ 



n 



n/2 



tr, 



e-^Q- / dnW^il - n)a^^''W,in) 



M 



i T2 

+ I dT2 /"rfriT(l-r2)e-(^2-"^)^T 



(n; 



where 



W,{r) 



X^e 



~tH 



T{t) = Ke 
Finally, by using the Duhammel formula 



-tH 



1 

dA= I dse^^-'^^{dA)e 





sA 



(5.2^ 



(5.29) 
(5.30) 



(5.31) 
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we can compute the derivatives of the exponential e 



-tH 



d,e-^ 



H 



ds e-'^^-'^^H^e-'" , 



(5.32) 



where 



H^ = a''^,^Uf3. (5.33) 

We can use this formula to compute the second derivatives of e~'^^ needed for the local 
coefficient ai 



T 





T S2 

+ ds2 dsi 





g-(s2-si)i/^^g-si//^ -(r-S2)// 



.e-^r-S2)H ^^^-{s2-s^)H -^^^-s^H 



where 



By isolating the overall exponential factor we get 



df,e 



-tH 



-P,ir)e 



~tH 



(3^{t)= I dse-'^'H^e'" , 



where 



which can be presented in the algebraic form 

M-) = —r. — H, = y: 



Ad 



H 



k=l 



-^[H,...[H,H,] 



k-l 



Here Adn is the operator defined by 

AdffA = [H, A] . 
By using the above formulas we obtain 

W^,(r) = [C,-/?,(r)]e--^ 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



(5.3^ 



(5.39) 



(5.40) 
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-tH 



where 



T{T) = i[a + 2r(3,{T)]e 



a = c,r - re, - iX'r 



(5.41) 

(5.42) 
(5.43) 



Thus, finally, summarizing this section we formulate the results in form of a theo- 



rem. 



Theorem 1 The heat trace of the operator L (\3. 6(\) has the asymptotic expansion as 
TiL2 exp(-tL) = (47rt)-"/2 [^^ + ^^^ + o{t^^ 



(5.44) 



where 



^0 = I dx / ^^ tr y e ^ , 



vr 



n/2 



(5.45) 



A-/ 



Ai = I dx 



d^ 



tr^e 



-"> -Q 



M M" 

1 



- / dn[C^ - ^^(1 - ri)]a^^[C - /3.(ri) 



1 r2 



+ /"rfra /"rfrie("2-"^)^[a + 2^,(1 - r2)Jle- 



(r2-Ti)H 





x[a + 2J'^/?^(ri)]L 



(5.46) 



ifere C^, /3^, a anc? J'' are given by eqs. ( \3.64\) , ( \5.33i) , (^.3%) , (^-W and ( ^.jSj) . 

5.2 Commutative Limit 

Let us compute the coefficient Ai in the commutative limit. We let 

af^'^ = gf^'^l ^ p = g^/^eH, Q = ql, (5.47) 

where g'^'^ is a nonsingular matrix, g = {detg'^'^)~^ , and is a scalar function. The 
operator L (|3.6CI| ) has the form 

L = -g-^/\d, + B, + <^^;)g'l^g^\d, + B, - (^,.)g-''' + g 

= -g-"\d, + B,)g''^g^^^{d, + B,)g-''' + g, (5.48) 
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where 

g = g + n0 + (7^'^0>,,. (5.49) 

and 

n = g-^'^d^g"^g^^d,. (5.50) 

Then 

H = g^^^i,i,, H, = g-^,,Up, J^ = g^''i., (5.51) 

C^ = -0,M-^(log^),/. + ^M, C^, = -(p,^,--{\ogg)^^,-B^, (5.52) 

« = -i,9'\. - 2e^(7^^S. , a = -i^g^\. + 2e^^'^'^i3. , (5.53) 

i5^{T)=TH^ = Tg^\^i^i^. (5.54) 

The integrals over the cotangent vector ^ are simply Gaussian and can be easily 
computed. First of all, we have 

'^^ e-"" = g'l\ (5.55) 



TT 



n/2 



So, we immediately obtain the coefficient A^ as the Riemannian volume of the manifold 

Ao= \ dxg^l^ . (5.56) 



M 

Next, we introduce the notation for the Gaussian averages 



(/>=/:|5J-"'^-''/K). (5^57) 



TT'' 



Then the Gaussian average of an exponential function gives the generating function 

(exp(^a;)) = exp ( -^g^yX^x"" \ , (5.58) 

where g^y is the inverse matrix of the matrix g^^ . Expansion in the power series in x 
generates the Gaussian averages of polynomials 

(eM.---eM.„+.> = 0, (5.59) 

(2n)! 
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where the parenthesis denote complete symmetrization over all indices. In particular, 



(1) = 1, i^f^^u) = ::g^,u, 



(5.61) 



i^fi^uCa^ls) — -iQinuQaP) 



jigfiuga^ + gf^agu^ + gtipgua) , 



15 



{^fi^uC,a^f3C,pC,a) — -^g(p.ugapgpa) 

_ I, 

— '^{gpugafsgpa + gpagutsgpa + gp^guagpa + g^pgafsgua + gpagafsgpu 
+gpugapgi3a + giiagiypgi3a + gppgupgaa + gppgavgpu + gp.ugapgi3u 
+gpvgaag/3p + gpaguagpp + g^pguagap + g^pgaagpu + g^agaug^p) ■ 



We have 



1 
Ai = j dxg^/Hiyl-q - j dn[C, - (1 - n)H,]g^^-{C, - nH, 

M 

1 T2 

+ dT2 dTi[a + 2(1 - T2)h]{a + 2ri/i) 







where 



Evaluating the integrals over the parameters Ti and T2 we obtain 



Ai = dxg^^^tiv(-q 

M 



-,pv 



2J^'r)B„B, 



pi^v 



-g'^'P^p'P,. - g'"" 



Hp + -{\o^ 



-i(7^^iJ^( log (?),,- ^^^'^( log (?),^( log. 



--g^^'H^H^ + ^A^ - -h\ + -/iM , 
6^ ^ 2 3 6 ^' 



(5.62) 



(5.63) 



(5.64) 



(5.65) 



(5.66) 
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where 

A = U''',- (5-67) 

Lastly, we need to evaluate the Gaussian averages. By taking into account the 

formulas above we obtain 
1 






flU 



^a/3^"^,M 



^logi 



,M ' 



2- 



(hX) 



g,.g'\ag''^,p, 



^-^{\ogg)^pg^\^^]^gp,g^\^g^\^, 



15 



^^> = ^9{p^vgc.^gpa)g''^g g^\^g 



-fp„St7 iiu a(3 



1 



(7^^(log(7)„(log(?),5-4(?^^^(log(7),. 



+2g^^g,pg'^'g>^\,g'^P^s + %a^^^.^"^,/3 + %.^'^^„^"''„ 
Here we used the equation 

(log^),^ = -gaiig'"^,p = g'"^ga(3,p 

We see that Ai does not depend on B at all and we get finally 



Ai = I dxg^'^N\-q-g^\ 

M 



,HY,y 



+ ^9"'{ ^ogg)A log (?),. + ^( \ogg),,g''\. 



(5.68) 
(5.69) 
(5.70) 
(5.71) 
(5.72) 

(5.73) 



(5.74) 



+-g(^ug'\ag''^,. 



1 

24 



g^agu/sg g ,'yg ,s 



(5.75) 



One can further integrate by parts here but we will stop here. Since the part that 
depends only on g must be an invariant, it must be an integral of a scalar quantity, 
which can be only the scalar curvature (there is only one such scalar). It turns out that 



Ai= dxg^'^N 

M 



-1 - S^'fef ,. 



6 



M 



(5.76) 
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where R is the scalar curvature of the metric g and q is defined by (|5.49| ). This coincides 



with the standard results [26, IS, 21] for the Laplace operator. 



6. Noncommutative Gauged Gravity 

Our main idea is that gravity has some new degrees of freedom. Instead of describing 
gravity by one tensor field we propose to describe it by new dynamical variables: the 
matrix valued tensor field a'^" and the density matrix field p. 

6.1 Matter in Gravitational Field 

The motion of a massive particle in the gravitational field should be determined by the 
principle of the extremal action. In general relativity the action is simply the interval 
S = —m J \/ —ds^, which means that particles propagate along the geodesies of the 
metric g^y. The parameter m is the mass of the particle. In our modified theory, we 
assume that a particle "splits" into several parts which propagate separately and have 
different masses. In other words, we propose to describe the motion of a particle in 
gravitational field by the action 






i=l 



'S'particle = - /^ "^(i) / J ^^^fi) ' (^-l) 



where the intervals dsf--. are defined by (|2.37 ). In the commutative limit the sum of the 



mass parameters determines the usual mass 

s 

5^m(,). (6.2) 



m 

i=l 



One can develop the whole Hamilton- Jacobi formalism starting from this for the particle 
mechanics based on Finsler geometry. We will not do it here but refer the reader to 
Rund's book |^ (in particular, Chap. 7 of the Russian edition). 



The dynamics of matter fields in gravitational field can be now described by the 
differential calculus developed in Sect. 3. A typical action for the matter fields has the 
form 

5ficid = [dx{- {ip, Lip) + W{ip)} , (6.3) 

M 

where L is an appropriate second-order differential operator constructed by the methods 
of Sect. 3 and W is a. potential. The first-order operators of Dirac type needed for 
spinor fields can be constructed similarly. 
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6.2 Equations of Gravitational Field 

Our goal is to construct an action functional that is invariant under both the diffeo- 
morphisms and the local gauge transformations. We also require that it is a noncom- 
niutative deformation of the Einstein general relativity, which means that our action 
should reduce to the standard Einstein-Hilbert functional in the commutative limit. 

We construct the classical action following the ideology of "induced gravity" as the 
large mass limit of the effective action. The classical action of gravity is then identified 
with the first two coefficients of the asymptotic expansion of the effective action as m — > 



oo (|5.8| ). In other words, we use the coefficients Aq and Ai to construct the invariant 
action functional of the noncommutative gravity. For a dynamical theory we need an 
invariant action functional that depends on the first derivatives of dynamical degrees 
of freedom. Such an invariant is given by the heat kernel coefficient Ai. Therefore, the 
invariant action functional of noncommutative gravity is a linear combination of the 
first two heat kernel coefficients 

^ = TTTT^^ |6 ^1 - 2A Ao] + ^^atter , (6.4) 

iovrG iV 

where ^matter is the action of matter fields G and A are the phenomenological coupling 
constants (Newton constant and the cosmological constant), and the coefficients are 
chosen in such a way that it has the correct commutative limit. This action is in 
some sort unique and the dynamical model described by it can be called the Induced 
Noncommutative Gauged Gravity. 

It is worth indicating the general form of the action. It can be represented sym- 
bolically as follows 

'5'= ^j,T / dxtTv\F{a,p)dada + F{a,p)dadp + F{a,p)dpdp 
IGttGN J I 

AI 

+F{a){6Q + 2A) + F(a, p)dB + F{a, p)Bb} + Scatter , (6.5) 

where d denotes the derivatives and F{a, p) denotes the coefficients that can only 
depend on a and p. Of course, all the coefficients are different and one has to include 
similar terms with all possible orderings of noncommutative factors. 
The fundamental equations of this model are 

Notice that the gravitational action does not depend on the derivatives of the vector 
field B^ (after integration by parts). Therefore, if the matter action S'matter does not 
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include a kinetic term for the fields B, then the variation with respect to B gives just 
a constraint, which simply expresses B in terms of the derivatives of the matrices a'^^ 
and p. Also, the matrix Q is fixed as a given self-adjoint positive definite function of 
a and p and their derivatives. Its form can be adjusted by imposing some additional 
physical requirements. This could be important for spontaneous symmetry breakdown. 
For simplicity, it can be just set to zero Q = 0. 

If we introduce a deformation parameter x according to ( 5.47 ): a'^'^ = gf^'^I + x/i^'", 
p = 5f^'^exp(I0 + >cip), and Q = ql + >cP, where h^'^ is a matrix- valued tensor field, 
and q are scalar fields, and ip and P are matrix-valued scalar fields, then in the limit 
X — >^ the action becomes 

^(0) = Y^ f dx g'/^ {R-2A- Gg^''^,^,^,^ - 6g) + 5^atter(o) , (6.7) 

M 

which describes the Einstein general relativity and a scalar field. Note that in the case 
q = —g^'^(t),^<p,y — □(/> the scalar field disappears. Other interesting choices of q are 
q = m^ or q = A(0^ — -m^)^. 

It would be certainly very interesting to compute the deformation corrections to 
the action ( p.7|) , that is to understand the precise way under which the commutative 
action ( |6.7] ) gets deformed to the non-commutative one (|6.4|). In general, we should 
obtain 

oo 

S{a, p, B) = ^(0) [g, 0) + J^ x'S^u) {g, <P, h, ^, B) . (6.8) 

fc=i 
Notice that the zero-order action S'(o) does not depend on the fields B. By a simple 
rescaling —>■ xcj) and q —* xq we could easily shift the dependence on these fields to 
the higher orders as well so that the zero order term 5'(o) [g) is just the Einstein- Hilbert 
action. The coefficients S^k) are polynomial in the fields 0, h, B and ip- Their general 
form can be read off from (|6.5| ). In principle, one could get the coefficients S(_k) from 
the general result of Theorem 1. However, such a calculation (even if straightforward) 
presents a real technical challenge and would require a separate paper. We plan to 
carry out this calculation in the near future and present the results elsewhere. 

7. Discussion 

In conclusion we list some interesting open problems in the proposed model. 

Uniqueness of the noncommutative deformation. First of all, it is very important 
to understand whether the proposed deformation of the general relativity is unique. If 
it is not unique, then what additional physical conditions should one impose to make 
such a deformation unique. 
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Interaction of noncommutative gravity with matter. One needs to find a consistent 
way to describe the interaction of the ordinary matter with gravity. In particular, 
to understand which physical matter fields should interact with the noncommutative 
(gravicolor) degrees of freedom and which should only feel the graviwhite part. 

Classical solutions and singularities. It would be very interesting problem to study 
simple solutions of this model, say a static spherically symmetric solution, which would 
describe a "noncommutative black hole" as well as a time- dependent homogeneous 
solution, which would describe the "noncommutative cosmology". One should study 
the problem of singularities of classical solutions. If this model is free from singularities, 
it would be a very significant argument in favor of it. 

Spontaneous breakdown of symmetry. One needs to understand whether it is pos- 
sible to introduce the spontaneous breakdown of the gauge symmetry, so that in the 
broken phase in the vacuum there is just one tensor field, which is identified with the 
metric of the space-time. All other tensor fields must have zero vacuum expectation 
values. In the unbroken phase there will not be a metric at all in the usual sense since 
there is no preferred tensor field with non-zero vacuum expectation value. 

Quantization and renormalization. Our model is, in fact, nothing but a generalized 
sigma model. So, the problems in quantization of this model are the same as in the 
quantization of the sigma model. 

Semi- classical (one-loop) approximation and heat kernel asymptotics. We point 
out that the study of the one-loop approximation requires new calculational methods 
since the partial differential operators involved are not of the so-called Laplace type 
(nonscalar leading symbol). For example, even the heat kernel coefficients A2 needed 
for the renormalization in four dimensions is not known in general. 

High-energy behavior. We expect that the behavior of our model at higher energies 
should be radically different from the Einstein gravity since there is no preferred metric 
in the unbroken phase, when the new gauge symmetry is intact. 

Dark energy. It would be very interesting to study the question whether the new 
noncommutative degrees of freedom of gravity could be accounted for the dark energy 
in cosmology. 

Low energy behavior and confinement. One could expect the gauge {gravicolor) 
degrees of freedom to be confined within some short characteristic scales (say, Planck 
scale), so that only the invariants {graviwhite states) are visible at large distances. Then 
the metric and the curvature would be only effective characterictics of the spacetime 
at large distances. 

Noncommutative deformation of Riemannian geometry. The Einstein spaces are 
manifolds with the metric satisfying the vacuum Einstein equations (with cosmological 
constant), i.e. R^^ = h^g^j^p with some constant A. In other words, the Einstein metrics 
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are the extremals of the Einstein-Hilbert functionaL The study of Einstein spaces is a 
very important subject in Riemannian geometry. It would be very interesting to study 
the noncommutative deformations of Einstein spaces defined as the extremals of the 
functional of matrix gravity. This would also have deep connections to "noncommuta- 
tive spectral geometry" . 
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